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Q ■ Abstract 
O 

Using a collective potential derived on the basis of the Generator Coordi- 
nate Method with Skyrme interactions we obtain values for the compression 



modulus of 40 Ca which are in good agreement with a recently obtained ex- 
perimental value. Calculated values for the compression modulus for 16 are 
also given. The procedure involved in the derivation of the collective potential 
is briefly reviewed and discussed. 
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I. INTRODUCTION. 



Among the collective vibrational modes of nuclei, the isoscalar giant monopole resonance 
(GMR), also called breathing mode, has received considerable attention due to the fact that 
its energy is related to the nuclear compression modulus. This property corresponds, in the 
limit of an infinite number of nucleons, to the nuclear matter incompressibility, which is of 
great interest in areas such as heavy ion collisions and nuclear astrophysics. The states cor- 
responding to this type of collective motion have long been experimentally identified JIJ and 
various microscopic theoretical approaches have been developed in order to account for and 
improve the deduced values of the compression modulus In recent years, more accurate 

experimental determination of the GMR energies |J led to renewed interest in the determi- 
nation of the nuclear matter incompressibility. Widely ranged theoretical approaches have 
been employed to obtain the compression modulus of finite nuclei and of nuclear matter, in- 
cluding microscopic calculations with the Gogny interaction and with generalized Skyrme 
forces H, self-consistent Hartree-Fock plus Random Phase Approximation treatments with 



Skyrme interactions 0, relativistic mean field calculations |10[ , a Thomas- Fermi [jTTJ and a 
Fermi liquid drop calculations [|12[ . 



Motivated by the recent experimental results 0, and by the discussion of the compression 
modulus in finite nuclei presented in flI3| , we report on a microscopic calculation of the 
compression modulus for 40 Ca and 16 based on the Generator Coordinate Method (GCM) 
|l4|Hl6|. This method allows for the microscopic determination of a collective potential 
V(q), q being the appropriate collective coordinate for the breathing mode, from which the 
compression modulus can be calculated directly as 

* q=qo 

where A is the nuclear mass number and go corresponds to the minimum of V(q). This 
approach has been developed in its general aspects many years ago |17j and can be directly 
used in the present context. We show in continuation that the collective potential does not 
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coincide with the diagonal part of the GCM energy kernel due to non-negligible off-diagonal 
contributions; moreover, the correlations embodied in the choice of the generator coordi- 
nate introduce anharmonicities in the collective potential, with the result that the deduced 
compression modulus comes out higher than that obtained from the standard variational 
approach. 

In our calculations we have used Skyrme interactions without Coulomb forces. Values 
for the compression modulus were obtained for 40 Ca and for 16 0. In the case of the heavier 
nucleus there is good agreement with the available experimental data especially when the 
Skyrme interactions SII and SIII are used. The GCM results also give, in the case of the 
lighter nucleus, values for the compression modulus which are only moderately lower than 
the values obtained for the heavier nucleus. Skyrme forces without density dependence such 
as SV lead to values for the compression modulus which are appreciably lower in both cases. 

In section II the GCM and its connection with collective Hamiltonians is briefly reviewed. 
In section III the collective potential for the breathing mode and the microscopically de- 
duced expression for the compression modulus are derived. Values of the incompressibility 
are calculated for five sets of Skyrme parameters. Finally, in section IV we present our 
conclusions. 



II. GCM KERNELS AND COLLECTIVE HAMILTONIANS. 

The well known Generator Coordinate Method is based on an ansatz for the collective 
wave function which is set up as 

1^ (0) = J f ( a ) ® a) da = J f (a) \a) da, 

where the {$ (f, a)} (or, for short, the {|a)}) constitute a set of auxiliary nuclear many-body 
wave functions parametrized by a (the generator coordinate), related to the nuclear degree of 
freedom one wishes to describe, / (a) is a weight function to be determined variationally and 
f is a shorthand notation for all the nucleon spatial coordinates. One way of constructing the 
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<3> (f , a) is through the use of Slater determinants of the lowest single particle eigenfunctions 
ip (r; a) of an ad hoc a-dependent auxiliary potential. The choice of the a-dependence of the 
auxiliary potential should therefore suit the particular collective degree of freedom under 
consideration (e.g. an overall scale parameter in the case of the breathing mode). 

For a given nucleon-nucleon interaction, we can write the nuclear many-body Hamilto- 
nian H and use the Rayleigh-Ritz variational principle to obtain the so-called Griffin- Wheeler 
equation ]l5f for the weight funtion f(a) 



J [(a\H\a'} - E(a\a')\ f (a) da' = 0, (1) 

where E is introduced as a Lagrange multiplier to account for wave function normalization. 
In this equation the energy kernel (a\H\a r ) and the overlap kernel (a\a') are functions of a 
and a' only, the nucleon coordinates having been integrated out. 

In principle, the solutions of the integral eigenvalue equation ([!]) give a variational ground 
state and a set of collective excited states expressed in terms of weight functions fi(a) and 
the corresponding eigenvalues Ei. Technical problems related to the overlap properties of the 
auxiliary many-body functions \a) in many cases prevent the implementation of such a direct 
approach, however. A discussion of these difficulties and of general methods to circunvent 



them can be found in flq , |19|j , where it is shown that by applying general transformations to 
the Griffin- Wheeler equation the determination of the collective spectrum can be reduced 
to a standard eigenvalue problem. 



A convenient choice for the nucleon-nucleon interaction is the Skyrme force | 2~0fl , which 
contains two body and three body terms. The two-body term is 
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vi,2 (n, r 2 ) = t (1 + x P a ) 5 (fi - r 2 ) 

S (r*i —r^k + k' 5 (r*i — f^) 
+t 2 k' ■ 5 (fx - f 2 ) k 
+iW {B\ + a 2 ) -k' x 5 (n - f 2 ) k, 

where to, xq, h, t 2 and Wq are adjusted parameters, and P a is the spin exchange operator. 



The three-body part of the Skyrme force is simply taken as a zero range Wigner force 



v 1,2,3 (n, r 2 , f 3 ) = t 3 5 (fi - r 2 ) 5 (f 2 - r 3 ) , 

where £ 3 is an additional parameter. For even nuclei this contribution can be rewritten as 
the density dependent two-body interaction 

v 1,2,3 (n, r 2 ) = ^ (1 + P a ) 5 (n - r 2 ) p |^!l±Ia^ . 

Since we are interested in the monopole giant resonance spectrum of spherical nuclei, we 
choose an oscillator potential with the oscillator parameter (3 as generator coordinate, which 
is thus directly related to the nuclear radius. The auxiliary many-body functions \(3) will 
therefore be Slater determinants built from single particle oscillator states {ip\ (f; (3)} with 
oscillator parameter (3. The corresponding overlap kernel then reads [P],|J 

where T is 6, 36 and 120 for 4 He, 16 and 40 Ca respectively. In order to handle this overlap 
kernel it is convenient to perform a change of variable by writing f3 = /9oexp (a), where (3q 
defines a reference length scale and a is a new generator coordinate. With this choice the 
overlap kernel appears as 

(a | a') = sech T (a — a') 
which, for the heavier nuclei, can be approximated by the Gaussian 

(a\a') ~ exp 
The energy kernel can be put in the form 



T , a2 

— (a — a ) 
2 v ; 



(2) 



(a\H\a') = (a\a') x [Cisech (a — a') exp (a + a')] 
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+C 2 cosh 3/2 (a — a') exp | — (a + a' 
+C 3 cosh 5//2 (a — a') exp - (a + a') 
+C4 cosh 3 (a — a') exp [3 (a + a')] , 



(3) 



where the coefficients Cj, i = 1 to 4 are given by 



^ ^ 2 ^ _ 3 t 



_ _ 3ti + 5t 2 9ti - 5t 2 ~ _ „ h 

8 (2tt) 16 (2vr) (tta/3 ) 



and ai, 02, 03, 04, and 05 are constants resulting from the integrations over the densities, 
having therefore values specific to each nucleus. 

Attempts to solve directly the resulting Griffin- Wheeler equation run into instability 
problems in the determination of the weight functions / (a). They can however be circum- 



vented by means of an analytical procedure discussed by Piza and Passos [18]. It consists 



in transforming the Griffin- Wheeler equation by introducing the new set of states 

l*> = J 7\WJk) 

where U k (a) and A (k) are defined through the diagonalization process for the overlap kernel 

U\ (a) (a\a')U k > (a') dada =A(k)5(k- k') . (4) 



Since the states \k) are orthonormal and complete |18j we are led to the new equation 



/ [H (k, k') -E]g [k') dk' = 0, 



where the new energy kernel is written as 



H(k,k') = (k\H\k') = J ^^(a\H\a')^^dada', (5) 



and 

g (k) = A 1/2 (k) [ Ul(a)f(a)da. 



For translationally invariant kernels (i.e., depending on a — a' only) such as (@), the diagonal- 
ization d|]) is performed simply by a Fourier transform, so that A (k) is also a Gaussian, and 
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there are no longer harmfull divergences. The transformations can be directly performed on 
EQ-dl) giving the new energy kernel in terms of the Fourier variable k. 

At this stage, numerical calculations can be performed and the energies and wave func- 
tions of the breathing mode can be obtained P|,[nj]. The resulting collective spectrum, 
obtained using the Skyrme interaction without the Coulomb interaction, is known PJi7[] 
and needs no further discussion. 



III. COLLECTIVE POTENTIAL AND COMPRESSION MODULUS. 

The point we want to stress in this note is that a value for the nuclear compression 
modulus can be deduced from the energy kernel given by Eq.((^). Although at first sight one 
is tempted to associate a collective potential to the diagonal part of the GCM energy kernel 
(a\H\a'), we will show that off-diagonal elements also contribute significantly to the properly 
defined collective potential. Since the energy kernel has the overlap kernel (a\a') as a global 
factor (see Eq.(|3])) and the overlap kernel is narrower for heavier nuclei, the contributions 
due to off-diagonal elements will be relatively more significant for lighter nuclei. 

In order to extract a collective potential from the transformed energy kernel H (k, k'), Eq. 



(||), we follow a method presented many years ago |I7[]. The first step consists of performing 
a double Fourier transform on H (k,k') 

H(k,k') -+H(x,x'). 

The resulting non-local energy function is then subjected to a Weyl-Wigner transformation 

h (q,p) = J(q - ^\ H \q + |) exp (iPijA dcr, 

where we have introduced the new variables 

x + x' , 
q = , a = x -x. 

It is then clear that the nonlocality of H (x, x') gives rise to the momentum dependence of 
the collective Hamiltonian. If we expand in the nonlocality parameter, a, we obtain the 
series 



O , ,,, O . 



# (?, cr) = (q--\H\q+-) = Y: # W (?) 5 H (a) , 

Z Z n=0 

where the coefficients are n-th moments of the energy kernel 

H (n) ( ) = L_l_ / # (g ff ) ^ 

n! j 

The resulting Weyl-Wigner energy function 

h (q,p) = £ f 5 (n) (a) exp (i P j) da 
n=0 "* \ nJ 

is associated to the collective Hamiltonian operator 

OO ( _A\ n 

H (q,p) = H&> (q) + £ {{.•■ W (?) ,p}} -} 

in which the p-independent term (q) is interpreted as the collective potential, namely 

V (q) = # (0) {q)= J H (q, a') da'. 
We will not discuss the inertia associated to the breathing mode, since this has been already 



before [17]. In what follows, we will concentrate instead on the collective potential. 

In order to evaluate V (q) let us return to the GCM expressions and expand the reduced 
kernel H (a, a') /N (a, a') (which can be read directly from Eq. (|3|)) as follows. Introducing 
the variables rj = a — a', and 7 = ^4^-, the general term of the series expansion of this 
object around the minimum 70 of its the diagonal part is 

H (a, a') 



1 <9 nH 

CWT (7 - TO)" 



N(a,a') 



- 1 7=70 ,v=o 



n\m\ dri m dY' 

As a result of the symmetry of the GCM kernels only even values of m and n appear in this 
expansion. The collective potential V(q) is correspondingly also given as a sum of terms 
Gm n (q — 70) which can be expressed as 

G mn (q ~ 7o) = ^ / exp rj-d V J J ^^"^ (7 - 7o)^7- (6) 

This is a somewhat symbolic expression, the meaning of which is specified by the prescription 
that the 7 integration is to be performed first in terms of derivatives of the delta function 



S(k), so that the integral over k gives derivatives of the /c-dependent part of the remaining 
integrand evaluated at k — 0. Following this procedure one may collect the powers of (q — 70) 
and write the collective potential as 

V 

Numerical results indicate sufficient convergence when this sum is truncated at v = 6 in the 
two cases studied here. The relevant coefficients V u are given by 

Vq = Coo - -^rpC 2 + j^C 2 o + Yof2^ 04 + f2^ m ~ ~^P 22 _ 64T3^ 06 + j^P m 
+ 16T3 C ' 24 ~ ir^ 42 



V 2 — C 2 - 2fC 4 + -jPl2 + ^06 - 7ypC24 + 7^642 (7) 



5 1 
V 3 = C 3 - ^^05 + C23 



15 1 
V 4 = Coa - ^Co6 + j, C24 



05 



Ve; — Co6- 

The above results show explicitely that the collective potential does in fact contain 
contributions from the off-diagonal terms of the GCM energy kernel, as mentioned before. 
In fact, the diagonal part generates only the terms Co n . Furthermore, it is also evident that 
the minimum of the collective potential does not coincide with that of the diagonal part 
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of the GCM energy kernel due to the off-diagonal contributions, which therefore introduce 
corrections to the equilibrium radius of the nucleus. Furthermore, we obtain the main result 
that the compression modulus, defined by 

1 d 2 



K 



Adq 



2 V(q) 



Q=Qo 



where go is the minimum of the colletive potential, does not coincide with the value obtained 
when the standard variational principle (based on the use of the diagonal part only of the 
GCM energy kernel) when the oscillator parameter is taken as the variational parameter. 
In fact, we obtain from the colletive potential 

2 r 



K 



A 



V 2 + W 3 (q - 7o ) + m (q - loY 



g=qo 



+ iov 5 (g-7o) 3 + 15^ (g-To) 4 " 

whereas on the basis of the simpler variational method we would have 



(8) 



(9) 



Comparison with Eq. (^) reveals the off-diagonal contributions to the predicted value of the 
compression modulus. 

Numerical results for 40 Ca using Eq.(|]) are shown in Table I using five different sets 
of Skyrme parameters, labelled as SI to SV. The values corresponding to the standard 
variational procedure (Q) are also shown for comparison. 

Table I - Compression modulus for 40Ca 



Compression \Interaction 


SI 


SII 


SIII 


SIV 


SV 


K (MeV) 


253.84 


229.08 


239.04 


217.86 


205.37 


(MeV) 


245.87 


222.28 


231.73 


211.57 


199.64 



It can be seen that the compression modulus obtained from the collective potential is higher 
than that calculated through the simple variational method. In fact, it is interesting to 
observe that the results obtained with interactions SII and SIII are in good agreement with 



the experimental value presented by Youngblood Jl3| 
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Using the sixth-order approximation to the collective potential we can also calculate the 
compression modulus for 16 0. The results are shown in Table II below. 

Table II - Compression modulus for 16 Q 



Compression \Interaction 


SI 


SII 


SHI 


SIV 


SV 


K (MeV) 


242.02 


213.78 


223.93 


202.60 


190.48 


K v (MeV) 


216.15 


192.17 


200.56 


182.71 


172.47 



In this case the value of the compression modulus obtained from the collective potential is 
about 10% higher than that obtained from the simple variational calculation. This indicates 
a theoretical value for the compression modulus of 16 which is not much lower than the 



observed experimental values for heavier nuclei [yj. It is also worth noticing that for 
SV, which does not include the three-body force and therefore does not present a density 
dependent term, the compressibility modulus is lower in both cases. 



IV. CONCLUDING REMARKS. 

Values for the compression moduli of finite nuclei have been obtained from a collective 
potential associated to a description of the breathing modes based on the Generator Coor- 
dinate Method. The procedure involved in the derivation of the collective potential from 
the GCM kernels has been reviewed. It has been shown that this procedure leads to a 
collective potential that embodies contributions coming from off-diagonal elements of the 
GCM kernels. They are associated to correlations introduced by the choice of generator 
coordinate and by the use of the Griffin- Wheeler ansatz, which lead to results which differ 
from those obtained from simpler variational procedures. In particular, values of nuclear 
radii and nuclear densities are different in both cases |22 . 



This approach has been applied to the calculation of the nuclear incompressibility for 
40 Ca and 16 0, using Skyrme effective interactions without Coulomb forces. The results 
are, as expected, higher than those obtained from simple variational approaches and, for 
interactions SII and SIII, are in good agreement with the experimental results given by 
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Youngblood et al |13| for 40 Ca. Our results for the lighter nucleus 16 O give values which are 



not too low when compared to the available experimental values for heavier nuclei and to 
the mean value of 231 ± 5 MeV assigned to nuclear matter by Youngblood. 

Acknowledgment: D.G. is partially supported by Conselho Nacional de Desenvolvimento 
Cientffico e Tecnologico (CNPq), Brazil. 
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